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Abstract 

A new method is proposed to calculate wave functions in /c-r-factorization in [1] as a comment 
about our paper [2] . We point out that the results obtained with the method are in conflict with 
the translation invariance and depend on the chosen contours for loop-integrals. Therefore, the 
method is in principle unacceptable and the results with the method cannot be correct. 

Recently we have pointed out in [2] that the /cr-factorization performed with off-shell partons 
is in general gauge variant and violated, because the perturbative coefficient function H contains 
gauge-dependent light-cone singularities and gauge-dependent ^-dependence. The /cy-factorization 
for the transition 7r° + 7* — > 7 has been examined in [3] in Feynman gauge at one-loop level. It is 
claimed in [3] that the perturbative coefficient function H calculated with off-shell partons is gauge- 
invariant. A proof of the gauge invariance at any-loop level is given in [3], where one shows the 
gauge invariance at N + 1-loop level with the induction method by assuming the gauge invariance 
at iV-loop level. But, the gauge invariance of the one-loop is not shown. This motivates us to 
study the problem in [2]. 

To determine the perturbative coefficient function H in the ^-factorization one replaces tt° 
with an off-shell quark pair. With the quark pair one calculates the form factor and the wave 
function. H receives contributions from the form factor and the wave function. In [2] it is shown 
with covariance arguments that the form factor does not contain the gauge-dependent light-cone 
singularities. The singularities come from wave functions. In the updated comment made by Li 
and MishimapQ a new method is proposed to calculate wave functions with off-shell partons. It is 
found that the gauge-dependent light-cone singularities disappear. In this updated reply we point 
out unacceptable flaws of the method. 

The definition of the wave function for 7r° moving in the z-direction is [2j [3] : 

<j>(x,k T ,C,») = J ^^e- p+2 ~-*^ (1) 

We have taken the light-cone coordinate system in which 7r° has the large momentum component 
P+. L u is a gauge link along the direction u^" — (u~^,u ,0,0) with u and C, is defined as 

C, 2 ~ A{u ■ P) 2 /u 2 . In the definition there is the momentum conservation in the +- and transverse 
direction due to the translation invariance. Because of this and z + = all — components of 
momenta in loops are integrated over from —00 to 00. It is noted that <j) only depends on x,kx,C 
and \i as indicated in the above. 

To extract the one-loop coefficient H^- 1 ' one needs to calculate (f> of an off-shell parton pair instead 
of 7T° at one-loop, as required in the fcj-factorization[3J. Then one calculates the convolution of 



with #(°), i.e., ® as a part of contributions to The mentioned singularity in 

comes from Fig. 2a, 2b and 2c in OH], where one has an integral of q~ from — oo to oo. q~ is the 
component of the momentum q carried by the exchanged gluon in the figures, the component q + 
and q± are fixed in <f> by the momentum conservation. The integral can be performed by a standard 
contour integral in the q~ complex-plane. One can takes a closed contour like those indicated in 
Fig. 5 of pQ for —R < q~ < R at first step by adding semicircles. After the integration one takes the 
limit R — > oo. In this limit one finds that all contributions to (j) from semicircles in Fig. 5 of pQ are 
zero. The remaining contributions to (j) are only from poles in the complex plane and contain the 
light-cone singularity [21 1]. In the limit the contribution fob from Fig. 2b in (2j [1] to <f> is nonzero 
only for < q + < kf with = (kf, 0, fcu_) as the momentum carried by the off-shell quark. 

In PQ it is suggested that one should keep these contributions to <fi from the semicircles to 
calculate the convolution (f)^> ® H^°\ which contains integrations over x and kx- x is related to q + 
through xP + = kf — q + in [21 [T] . After finishing the integrations in the convolution one then takes 
the limit R — > oo. In this way the light-cone singularity disappears because the contributions from 
semicircles given in Fig. 5 PQ give nonzero contributions to the convolution in the limit. Since the 
limit is taken after the convolution, this implies that <j> depends on an extra parameter R which 
is in fact a cut-off for the large \q~\- The existence of a finite R actually is in conflict with the 
translational invariance. We take the contribution fob as an example to illustrate this in more 
detail in the following. 

It is well-known that the wave function in Eq.(l) is zero for x < and x > 1. The boundary 
x = 1 can be obtained by inserting a complete sum of intermediate states in Eq.(l) and using the 
translational invariance. With the finite R fob is not zero for q + < 0. This nonzero contribution 
is crucial to cancel the light-cone singularity of fob with q + > OpQ. If one takes the limit R — > oo 
before the convolution, fob is zero for q + < 0. From the translation invariance one can show that 
fob or contributions from similar class of diagrams are zero for q + < without any calculation. 
The contributions from Fig. 2b in [21 1] and its higher-order corrections, in which there is no gluon 
attached to the antiquark line and the gauge link L u (oo, z), are actually proportional to the contri- 
butions to a distribution of a single quark with the momentum ki . The distribution can be defined 
by replacing 7r° in Eq.(l) with the single off-shell quark and delete some operators in Eq.(l) as in 
the following: 

dz_^ etyktz -„ l?± r kT{ol ( L t (oO)0) _ q{zMh)) \ z+=Q . (2 ) 

fob is proportional to one of the leading order contributions of the above distribution with xP + = 
ykf. Inserting a complete sum of intermediate states and using the translation invariance, we have 
that the above distribution must be zero for y > 1. Here y is related to q + through yki = k~l — q + . 
One can conclude without any calculation that the above distribution and fob must be zero for 
q + < 0, or the contributions from similar diagrams are zero if the out going quark line carries 
a +- momentum component is larger than that of the incoming quark line. This conclusion is in 
agreement with the physical expectation that partons entering hard scattering have positive +- 
energy. Therefore, fob ^ with q + < is in conflict with the translation invariance. To avoid the 
conflict one has to take the limit R — > oo before the convolution. This results in fob = for q + < 0. 

In fact there is certain ambiguity in fob if one does not take the limit R — > oo. In the method, one 
actually changes the order of the integrations in the convolution fob^H^ with the q~ -integration 
in fob by implementing certain prescriptions, where one first deforms the contours along the real 
axis in the g~-plane for the ^"-integral in <p to those in Fig. 5 of pQ. With these contours one finds 
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the nonexistence of the light-cone singularity. If one flips the closed contour in Fig. 5b of pQ, i.e., 
by changing the contour in the upper-half q~-plane into the lower-half q~-plane, one will find a 
different result from Fig. 5b than that given in Eq.(12) of [1] by noting that fact: The contribution 
from the single pole in the upper-half g~-plane is the same as that from the double pole in the 
lower-half (/"-plane, and with the flipped contour the corresponding integral of 8 from tt — > in 
Eq.(12) of [TJ is changed to that of 9 from —tt —* 0. The final results with the flipped contour will 
still contain the light-cone singularity. The cancelation of the light-cone singularity only happens 
with certain contours like those given in Fig. 5 of pQ. This is in conflict with the general expectation 
that physical results should not depend contours used for loop-integrals. 

It is interesting to note that one will not find the singularity in the convolution if one simply 
changes the order of the integrations of 4>2b <8> by ignoring the constraint of allowed t/ + -regions 
from the translational invariance, where one first integrates over q + by taking the region with 
q + < into account and then perform the integration over q~ . However, from the translational 
invariance we know that the region with q + < should be excluded in the integration of q + . Taking 
this into account, 4>2b and its convolution with has the mentioned light-cone singularity. 

The above discussion tells that the final results with the method will depend on chosen contours. 
The method implies a modification for the wave function of parton states through implementing 
a cut-off and specially chosen contours. It is unclear how to implement the modification to the 
wave function at hadron level, i.e., to that in Eq.(l). With the implement the calculated wave 
function is in fact no longer the wave function employed in the /^-factorization m Eq.(l). Since 
the results obtained with the method are in conflict with the translational invariance and depend 
on the chosen contours, the proposed method is in principle unacceptable. 

To conclude: In the suggested method of [1] an extra cut-off is introduced into wave functions 
and special contours for — components of loop momenta are taken. The results obtained with 
the method are in conflict with the translation invariance and depend on the chosen contours. 
Therefore, the method is in principle unacceptable and the obtained result with the method can 
not be correct. 
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